Gaussification and Entanglement Distillation of Continuous- Variable Systems: 

A Unifying Picture 
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Distillation of entanglement using only Gaussian operations is an important primitive in quantum commu- 
nication, quantum repeater architectures, and distributed quantum computing. Existing distillation protocols 
for continuous degrees of freedom are only known to converge to a Gaussian state when measurements yield 
precisely the vacuum outcome. In sharp contrast, non-Gaussian states can be deterministically converted into 
Gaussian states while preserving their second moments, albeit by usually reducing their degree of entangle- 
ment. In this work - based on a novel instance of a non-commutative central limit theorem - we introduce a 
picture general enough to encompass the known protocols leading to Gaussian states, and new classes of pro- 
tocols including multipartite distillation. This gives the experimental option of balancing the merits of success 
probability against entanglement produced. 

PACS numbers: 03.67.Bg,42.50.Ex 



Entangled quantum states are the fundamental resources 
that enable quantum key distribution, quantum communica- 
tion, and instances of distributed quantum computing. Real 
physical systems are affected by decoherence and non-ideal 
apparatus that degrades the quality of experimentally prepara- 
ble quantum states. However, entanglement distillation pro- 
tocols provide a means of converting many copies of par- 
tially entangled states into a smaller number of more entan- 
gled states [1J. When entanglement is required over very long 
distances, distillation can be implemented at regular intervals 
called repeater nodes (2J-H). Photonic systems that carry en- 
tanglement in continuous degrees of freedom are difficult to 
manipulate arbitrarily. However, so-called Gaussian opera- 
tions are more easily implemented by a combination of beam- 
splitters, phase shifters and squeezers. Furthermore, prepa- 
ration of Gaussian states is routine in many laboratories, and 
such states are especially useful for numerous quantum infor- 
mation tasks. 

Unfortunately, Gaussian operations are quite limited in 
their capacity to distill entanglement. In particular, a series 
of no-go theorems have shown that with only Gaussian re- 
sources, it is impossible to increase entanglement [5]. These 
results can be circumvented when non-Gaussian resources are 
available. Given an appropriate resource, Gaussian operations 
can simultaneously increase the entanglement and make the 
state more Gaussian, a process we refer to as the "Gaussi- 
fication protocol" (GP) |6| (for steps towards experimental 
realization of this and related protocols, see Refs. [71). Al- 
ternatively, the theorem can also be circumvented by using 
a non-Gaussian operation implemented by photon detectors. 
Entanglement distillation can be achieved by a combination 
of photon subtraction, a de-Gaussifying operation, and Gaus- 
sification HE]. 

While these techniques give hope for simple realizations of 
quantum information protocols, they do not exploit the rich- 
ness of Gaussian operations. Specifically, GP utilizes only 
projections onto the vacuum. These projections are feasible 
if reliable detectors are available that distinguish zero from 
one or more photons. However, strictly speaking these detec- 
tors do not fall within the realm of Gaussian devices. Per- 



FIG. 1: (i) A single step of the general class of protocols (GG) con- 
sidered, illustrated for three parties. This embodies the known Gaus- 
sifier (GP) entanglement distillation schemes based on projections 
onto the vacuum, or the extremality protocol (EP) mapping unknown 
states onto Gaussian ones with the same second moments. The Gaus- 
sian projection can be arbitrary (including (ii) vacuum projection, 
(iii) homodyning, (iv) eight-port homodyning, (v) tracing out); and 
schemes with finite widths of acceptance. 



forming eight-port homodyne detection [9| and postselecting 
on the vacuum outcome achieves the same projection, but 
this will have zero success probability when postselecting ex- 
actly on the vacuum measurement outcomes. In this work, 
we prove Gaussification for a wide class of truly Gaussian 
protocols with nonzero, and tunable, success probabilities for 
multimode states. All known feasible distillation protocols, 
including our protocols, are iterative and consume a num- 
ber of copies exponential in the iterations required. As such, 
our scheme's capability to increase success probability sig- 
nificantly improves the prospects of distillation and repeater 
implementations with only modest resources. Our techniques 
also open up the perspective of directly distilling into multi- 
partite Gaussian states. 

In addition to the practical applications of our results, our 
analysis provides a more intuitive explanation of the phenom- 
ena of Gaussification. In the original distillation protocols [6|, 
the process of Gaussification is quite mysterious, but it is very 
apparent in the protocol of Ref. iflOll which provides an al- 



ternative method that uses no measurements, referred to as 
the " extremality protocol" (EP), as it is used to show the ex- 
tremality of Gaussian states with respect to several proper- 
ties. Yet EP can convert many non-Gaussian states into a more 
Gaussian state while conserving the expectation value of ob- 
servables quadratic in position and momentum. Although the 
EP does Gaussify, its capacity for increasing entanglement is 
restrained by its deterministic nature. The proof of EP ele- 
gantly employs the central limit theorem [11] that explains the 
ubiquity of Gaussian distributions in classical statistics and 
nature itself. Our approach unifies GP and EP within a com- 
prehensive theory of "general Gaussification " (GG) founded 
on a non-commutative central limit theorem and so provides 
an intuitive mechanism for Gaussification. In addition to bi- 
partite entanglement distillation, our approach reveals whole 
new classes of protocols that we discuss. 

We consider GG protocols that can be implemented itera- 
tively, with the (n + 1) iteration as follows: 

1. Take two copies of an m-mode state p n shared between 
m parties; 

2. Each party applies a 50:50 beam-splitter transformation 
between their pair of modes; 

3. Every party makes a Gaussian measurement on the out- 
put of one beam-splitter port; 

4. The parties compare measurement results and postse- 
lect such that the operation implemented is Gaussian; 

5. The output state p n +i is used for the next iteration. 

Formal description of GG protocols. There are 2m modes 
involved in the protocol, and we label annihilation operators, 
cijk, with two indices; The index j = 1,2, ... ,m labels the 
respective party and k — 1,2 the copy at a particular node. 
The beam-splitters, in the Heisenberg picture, perform 



Ua j>k tf = (o,-,! + (-1)* 



*j,2 



(1) 



The measurements at step 3 can be homodyne, eight-port ho- 
modyne, or any other Gaussian measurements projecting onto 
a state II m for measurement outcome m. We are interested 
in Gaussian protocols that postselect on a set of measurement 
outcomes, and mix over all accepting outcomes, 



Pn+l OC 



j dmP(m)tr 2 [U(p n ® p n )U\l ® n m )] 



(2) 



where tr 2 denotes a partial trace over the second copy, with 
k = 2. The integral over measurement outcomes is weighted 
by P(m). The weights P(m) = and P(m) = 1 corre- 
spond to a rejection and an acceptance, respectively, but we 
also allow for probabilistic strategies where -P(m) gives the 
probability of acceptance. The protocol is described by a sin- 
gle operator we call the filter 



n = / rfmP(m)n r 



(3) 



We allow for arbitrary Gaussian filters, II, that are invertible 
and proportional to a fully separable Gaussian state with van- 
ishing first moments and finite energy. Certain interesting 



cases, such as GP, EP and protocols using precise homodyne 
detection are included as limits within this family of filters, so 
keeping full generality. With this notation 



Pn+1 oc tr 2 [U(p n <g> p n )U\t <g> II)]. 



(4) 



Phase space. Before we give our results, we review phase 
space representations for the position and momentum observ- 
ables of an m-mode system, labeled as 

R = (Pi; R-2, • • ■ ; P2m-1, P2m) = (-X"l > Pi, ■ ■ • , X m , P m ) , 

where Xj = (a'- + a,j)/^2 and Pj = i(al — Oj)/-\/2 for 
j = 1, . . . , m. The canonical commutation relations between 
the coordinates are embodied in the symplectic matrix S. The 
covariance matrix of an operator A records the second mo- 
ments of these observables. 

{T A )j,k = tr({Rj - {d A )j,R k - {d A ) k }+A), (5) 

where {., .}+ denotes the anticommutator and first moments 
are ((1a) j = tr(RjA). Furthermore, we make use of char- 
acteristic functions, \A ■ R 2m — >• C, that encode all the in- 
formation of A as Xa(t) — tr(D(r)A), where D(r) is the 
displacement operator, D(v) = exp(?r • R). Such a function 
is said to be Gaussian when 



Xa(y) = exp(ir • d^ - r T rAr/4). 



(6) 



If the operator A is a physical state, its covariance matrix will 
be real. However, an instrumental tool in our analysis is that 
we work with A = a n = p n II/tr(p„n). Indeed, we will not 
employ characteristic functions of states satisfying the con- 
ditions of Bochner's theorem ifTTI [121 . but of more general 
objects, hence leading to more general complex-valued func- 
tions. Since II is invertible, a n uniquely defines a quantum 
state p n . 

A new non-commutative quantum central limit theorem. 
With these definitions at hand we can state our first result, 
with a stronger form of convergence demonstrated later. 

Theorem 1 (Convergence of general Gaussifier protocols) 

Consider an initial state p = po, with associated operator 
a = pll/tr(pll) such that the following conditions are 
satisfied: (i) d a — 0; (ii) |xo-(r)| < 1 for all r; (Hi) the 
covariance matrix 

r Poo = (r n - *£)(r n - rvr^rn + *£) - r n> (7) 

exists and is positive definite. Let p^ denote the Gaussian 
state with covariance matrix r poo . GG with filter II causes p n 
to weakly converge to p^ in the following sense: If \x) and 
\y) are eigenvectors ofH, then for all e > there exists a n e 
such that for all n > n f 

\(x\p n \y) /tr(p n n) - (x|Ax>|y}/tr(pooII)| < e. (8) 

Proof of the statement. Much of the basic structure of the 
proof follows an argument of a non-commutative quantum 



central limit theorem for quantum states [10-13]. The prob- 
lem once one allows for Gaussian measurements is that the 
characteristic function of the output is unwieldy indeed. We 
circumvent this problem by a bold step: We put in an addi- 
tional filter II at the output "by hand", in order to exploit sym- 
metry, at the expense of having to consider \a n of different 
objects, <7 n , This will then lead to the desired result. From 
Eq. (pi we have, 

X«r„ +1 (r) <x tr[(D(r) ® l)U(p n <g> p„)*7 f (II <g> n)]. (9) 

Using the cyclicity of the trace and Eq. (fTh 

XCTn+1 (r) oc tr[D(r/V2)® 2 p® 2 U^I1® 2 U}. (10) 

Next we recall that Gaussian states with zero first moments 
commute with beam-splitters, [/IT® 2 = H® 2 U, such that 

X<r„ +1 (r) ex tr{[D(r/V2) Pn Tl}® 2 } ex tr[D(r/V2)a n ] 2 , 
= Xa Jr/V2) 2 = Xa (r/VN) N , (11) 

where in the last equality N = 2™ +1 , iterating the formula. 
The key to the simplicity of this formula is to consider conver- 
gence of a n rather than directly p n . By introducing an addi- 
tional projector within the trace, symmetry allows us to com- 
mute through the beam-splitter unitaries, which is the essen- 
tial simplifying step. To find the limiting characteristic func- 
tion, we consider a given phase space point r and the function 
fa„ +1 ,r ■ K ->■ C defined as fa n + lt v (t) = Xa„+i(^o)- In 
the spirit of a classical central limit theorem iflOl flTl [131 but 
for non-Hermitian operators we write 
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N 



,(12) 



which converges assuming that second moments are finite, 
/(0) = 1 and |/(t) | < 1 for all t. This last condition, 
which is always satisfied for classical characteristic functions, 
may be violated for non-Hermitian <r n . However, provided 
|Xo-(r)| < 1 f° r ai l r > we find ( see App. Ia} in the limit of 
large n that 



lim fa n+ i(t) ex exp( 



-t 2 r^T r 



ro/4), 



(13) 



pointwise in t. Setting t = 1 shows pointwise convergence of 
Xcr n+1 for each phase space point i"o. Furthermore, following 
the reasoning of Refs. lfT2l[T3l , this entails that any operator 
B that is absolutely integrable such that J \xB(?)\dr < oo, 
has a convergent expectation value ti(Ba n ) — > tr(_Bcr 00 ) (see 
App.lBJi. Setting B = \' y 1 \y){x\ where T\\y) = X v \y), then 
for large n 



(x\p n \y) (x\<7 n \y) {x\aoo\v) (x\p a 



tr(p„n) 



A 



a 



A 



a 



trip^U) 



(14) 



All that remains is to show that Gaussian a^ entails a Gaus- 
sian poo. First we observe (see App. O) that the product of 
two operators has a characteristic equation 



Xa x (q) 



X,,- 



3(i")xn(q-r)exp(-ir T crq/2)dr. (15) 




FIG. 2: The degree of entanglement of the initial state p in terms of 
the log-negativity versus the target Gaussian p^ for (i) bipartite state 
vectors |*a) and (ii) tripartite state vectors 1$^}. In both figures the 
variable < A < 1 controls the degree of postselection. The curve 
varying with A is the entanglement of p^ , and the curve constant in 
A is the initial entanglement of p. 



For Gaussian xn and x Poo , the integral is a multi-variate 
Gaussian integral that evaluates (see App.[C| to another Gaus- 
sian with the covariance matrix taking the form of a Schur 
complement ifTBl . 

r ff = r n - (r n + iS)(r Poc + r n )- 1 (r n - ffi), (16) 

Rearranging this formula for r Poc gives us the covariance ma- 
trix for the convergent state p^ in terms of T a and Tn as in 
the theorem. 

Examples. We have introduced a broad class of protocols 
for which our theorem indicates Gaussification. However, for 
concreteness it is helpful to keep in mind a simple class of pro- 
tocols. Consider when each party performs eight-port homo- 
dyne measurements that project onto a coherent state. When 
the fc th party obtains outcome a^ projecting onto coherent 
state vector \ctk), we declare the iteration a success with prob- 
ability P = exp(— \a\ 2 /2c 2 ). The degree of postselection is 
quantified by a real variance c. It follows that the filter is in- 
deed proportional to a Gaussian state with covariance matrix 
T n = A" 1 !, with A = (1 + c 2 ) 3 / 2 . This class of protocols 
is important as it contains GP and EP as limits A — ► 1 and 
A — » 0, respectively. 

We now consider the degree of entanglement that is 
achieved by applying our protocol with eight-port homodyne 
measurements. First we consider the well-studied bipartite 
state vector l 1 !^) ex |0,0) + A|l,l), and present the log- 
negativity [16 1 of the convergent Gaussian in Fig. 13] Vary- 
ing the parameter A interpolates between the entanglement 
achieved by GP and EP, with increased yield compensating for 
reductions in entanglement. We also analyzed a tripartite en- 
tangled state vector |$ M ) ex |0, 0, 0) + /j(|l, 1, 0) + |1, 0, 1) + 
|0, 1, 1)), having up to two photons in three modes, and for 
the log-negativity summed over all 3 bipartitions. In fact, in 
this way, one can straightforwardly engineer multi-partite hy- 
brid distillation protocols for quantum networks, giving rise 
to primitives in repeater architectures where entanglement is 
shared across many repeater nodes. These would overcome 
the known limitations providing a road-block against entirely 
Gaussian continuous-variable repeater networks II 1 71 . 

Strong convergence. Our previous theorem proves a con- 
vergence result identical to that of Ref. [6| and GP. However, 
it would often be preferable to have convergence of (x\p n \y) 
to (x\p OQ \y) in a stronger sense (without the addition factors 



of ti^pnll)" 1 and tr( j o oo n) _1 and as a convergence in trace- 
norm). For most physically relevant (see App. [D} instances 
of initial quantum states, we now show that this is indeed the 
case. We will make use of expectation values of normally or- 
der operators, 

c$* = tr{V(®" k l =1 a x k *)\®? =1 af)Vla n ], (17) 

where x, y £ N m and V is the Gaussian unitary such that 
yilV^ is a thermal state. 

Theorem 2 (Strong convergence) In addition to Theorem^ 
if for all x, y £ N m the expectations values of normally or- 
dered operators satisfy, |ctQ ,y | < a^y for a^y = \a^\. It 
follows that for all e > there exists an n e such that for all 
n > n e , we have \ \ p n — p^ | | < e, where \\.\\ is the trace norm. 

The conditions of the theorem are stated technically, but phys- 
ically prevent overpopulation of higher Fock numbers. E.g., it 
is easy to check these conditions are meet in all the low pho- 
ton examples analyzed in Fig [2] We begin by first showing 
tr(Ilp n ) converges to tr(IIp 00 ), which in terms of a is 



tr(rrV) = tr(rrVn)/tr(/9n) = tr (pn)' 



(18) 



The inverse filter, II -1 , has an exponential form that can be 
Taylor expanded and, using the bosonic commutation relation, 
normally ordered such that tr(n -1 <7„) = X) x v ^ x >y Q « Y ' 
where q Xyy > 0. Hence, when the conditions of our theo- 
rem hold, we conclude tr(n _1 cro) < ti^n^Voo), and we 
proceed by showing this holds for all a n . Iteratively we have 
(details in App. [EJ 



x.y 

*n+l — 



E 

u<x,v<y 



(19) 



where C x ,y . is a combinatorial quantity that is non-negative 
and real. Under our assumptions the absolute values obey 



*n+l 



< V C x ' y a u ' v a x ~ uy ~ 



» x.y 



(20) 



u<x,v<y 



so initially satisfying our conditions entails satisfaction for all 
n. Hence, for all n we deduce tr(n _1 (T„) < tr(n~ 1 cr 00 ) and 
equivalently tr(p„n) > tr(p oo n). Next we bound tr(p„n) 
from above. Consider a finite rank projector P that commutes 
with n, then B — U~ 1 P has an absolutely integrable char- 
acteristic function. Hence, for arbitrarily small 5 > there 
exists an rig such that for n > n$ 



tr(Ba n ) ~ tr(-Beroo) = 



tTJPpn) tl(P Poo ) 

tr(IIp n ) tr^poo) 



> S. (21) 



Since tr(Pp n ) < 1, and P can be chosen so ti^Pp^) — 
1 — 5' is arbitrarily close to unity, we conclude ti^Ilpn)^ 1 — 
tr(npoo) -1 > — e where e = 5 + 5'tr(p 00 H)~ 1 is again small. 
This inequality rearranges to tr(np„) < (tr(IIp 00 ) _1 — e) , 
giving an arbitrarily tight bound from above. Combined with 
our lower bound we conclude that tr(Hp n ), under the stated 



assumptions, converges to tr(n / o oc ). As such (x\p n \y) con- 
verges to (x\poo\y), and as is well known this entails trace 
norm convergence (See Ref. lfT2ll and ApplFjl. 

Summary. We have introduced a framework for construct- 
ing a range of new protocols for entanglement distillation and 
manipulation. At the same time, this work provides a unified 
framework for existing protocols leading to Gaussian states: 
Notably, the mysterious emergence of Gaussian states in dis- 
tillation schemes is once again related to an instance of a quan- 
tum central limit theorem, albeit for a much more broader 
class of protocols than previously considered. This framework 
also allows us to look at Gaussification in experimentally real- 
istic acceptance windows, and to trade-off different figures of 
merit against each other. Such trade-off control is essential as 
previous proposals are so heavily postselective that over sev- 
eral iterations the success probability would reduce dramati- 
cally, whereas our protocols offer an arbitrarily good chance 
success. Potential for future research is broad as a unique pro- 
tocol is defined by every separable Gaussian state. For ex- 
ample, our techniques can be applied to homodyne detection 
protocols for either multimode entanglement distillation, or 
single mode squeezing enhancement. 
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Appendix A: Central limit theorems reviewed 

Here we present additional details on the convergence of 
characteristic functions to a Gaussian function. Our proof 
follows the classical proofs, such as in Ref. [14|, but since 
cr = pn/tr(pn) is generally neither positive nor Hermitian 
there are some notable subtleties and technicalities involved. 
Ultimately, we aim to prove the following. 

Theorem 3 (General quantum central theorem) If f a ,r (*) 
satisfies |/o-,r o (0l — l/ or all t £ R then in the limit of large 

N 



WVVAO 



, N _v --"' / 4 



(Al) 



pointwise in t, with v as the appropriate second moment. 
Specifically, pointwise convergence means that for all t £ M. 
and all e > there exists a N e such that for all N > N e we 
have 



\U !T0 (t/VN) N -e- vt2 / 4 \< 



(A2) 



Setting t = 1 directly entails pointwise convergence of the 
characteristic function x„ n , and furthermore this is uniform 
on compact regions, such as a ball in phase space. 

Corollary 1 (Convergence in compact regions) Within any 
compact region R, \<y n converges uniformly to Xa^- That 
is, for any e > there exists an n e such that for all n > n E 
and all r £ R, we have \ Xa n ( r ) ~ Xo-oo ( r ) I — e • 



Pointwise convergence is uniform on compact regions when- 
ever the gradient is bounded within R for all n. To prove a 
gradient bound for all n, we first observe that for n = a 
bounded continuos function always has a gradient bound on 
compact sets. That is, an 77 can always be found such that for 
all r e R we have G (r) = ||Vx CT0 (r)|| 2 < rj\r\, where ||.|| 2 
denotes the 2-norm. This property is iteratively preserved for 
all a n , since 



G„(r) 



= ||V XCT0 (rMvn| 2 , 

= N\\[Vx, (r/VN)].[xa (T/VN) N - 1 ] 

< y/NG (r/VN)<ri\r\. 



Using the point in the compact region that gives the largest 
value |r| giving the desired gradient bound for all r e R and 
all n. This corollary will prove useful in App. [B] 

We prove several useful lemmas before directly addressing 
this proof. First we observe, 



U,r (t) = tv[D(tr )a], 

/OO 
dxe lxt F a , ro (x), 
OO 



(A3) 



where F aiTa (x) is the representation of a along a line in phase 
space defined by the unit vector r (e.g., the position or mo- 
mentum representation in case of a unit vector pointing along 
the axis of phase space). This is similar to a probability den- 
sity for measurement of a particular quadrature; however, be- 
cause a is neither positive nor Hermitian, the function F aro is 
not generally positive or real, which is the main cause for cau- 
tion. Yet it is sufficiently well behaved, with the two following 
lemmas proving useful. 

Lemma 1 (Convergence of absolute integrals) The abso- 
lute integral satisfies J_ x 2 \F a ^ {x)\dx < 00. 



Now clearly (n — ip)(H + ip) > 0, such that 

n 2 + p 2 > ipH - iHp 
and (II - p)(II - p) > 0, so that 



II 2 + p 2 > Up + pU. 



(A4) 



(A5) 



Similarly, II 2 +p 2 > -ipH + illpand H 2 + p 2 > -Tlp-pU. 
One finds after a few steps that, for any vector \ip),v/e have 

M P ii\<p)\ < |(v-|(p 2 + n 2 )|^)|/y2. (A6) 

Hence, using the triangle inequality, 

< i^wki^wi . (A7) 

V2|tr(pII)| 

That is to say, whenever the second moments of n 2 and p 2 are 
finite, the absolute integral converges to a finite value. 
The second lemma we require is as follows. 



Lemma 2 (Expansion of functions) For a given unit vector 
Yq and an operator a that is unit trace and zero displacement 
d CT = 0, the above function / CTiI . can be expressed as 



U, To (t)=l-\vt 2 K(t) 



(A8) 



where K(t) — > 1 as t — >• 0. 



This lemma asserts that higher order moments drop off suf- 
ficiently quickly for the purposes of establishing a central limit 
theorem. We make use of the following formulation of the ex- 
ponential function 



1 + itx- -t 2 x 2 L(x,t) 



where 



satisfying 



L(x,t) 



2 / (1 - u)e mtx du, 
'0 



\L(x,t)\<l 



(A9) 



(A10) 



(All) 



for i,(el. Furthermore, for any finite x interval, e lutx — > 1 
as t — >• uniformly in x. Consequently, we also have 
L(x, t) — > 1 as t — >• uniformly on any finite x-interval. Mul- 



tiplying Eq. i~A§\ by F aro and integrating w.r.t. x and using 
the unit trace of a and vanishing first moments we have 

1 f°° 

fa.ro (t)=l- 0^ / X2L ( X > V F °>*° ^ dX > (A12) 

^ J -co 



which we can rewrite as 



fa,r (t) = l--t 2 uK(t)dx, 

where v is the second moment defined by 



(A13) 



(A14) 



(A15) 



) 

and K{t) is the quantity satisfying 

/OO 
x 2 L(x 7 t)F^ T0 (x). 
-OO 

If we subtract Eq. < |A14| > from the above, we have 

/•OO 

v{K{t)-l) = 2 x 2 {L(x,t)-l)F a , ro {x). (A16) 



(A17) 



J — OO 

The absolute value satisfies 

/■OO 



/OO 
x 2 \(L(x,t)-l)F a!To (x)\. 
-00 

Noting that \L(x, t) — l\<2 since \L(x, t)\ < 1, and dividing 
the integral into two parts over some finite interval R and the 
complement R^- we deduce that 

u\K(t)-l\ < 2 [ x 2 \(L(x,t) - l)F airo (x)\ 
Jr 

+ 4 / x 2 \F^ ro {x)\. (A18) 



We now make use of Lem. ([T]i which tells us that the second 
integral always converges to a finite value. Therefore, for any 
e > the region R can be choose sufficiently large that the 
second integral is bounded such that 

v\K(t)-\\ < 2 [ x 2 \{L{x,t)~l)F a . T0 (x)\+^\9) 

JR 

Recall that for finite x interval, such as R, (L(x, t) — 1) uni- 
formly goes to zero with vanishing t. Hence, for any e > 
and x £ R there exists a t. such that for all \t\ < t e we have 



v\K{t) - 1| < 2e. 



(A20) 



We conclude that K{t) converges to 1 with vanishing t, prov- 
ing our lemma. 

We now employ the above lemmas to demonstrate point- 
wise convergence. We wish to bound the quantity 



6(t,N) = \f a , To (t/VN) N 
which can be expanded as 

vt" 



-ut 2 /A\ 



5(t,N) = (1 -'^Kit/y/N)) 



In the limit of large N we have that 

vt 2 — 



N 



,-^/i 



(A21) 



. (A22) 



5{t,N) 



1 



AN 



K(t/VN) 



N 



Vt 



AN 



(A23) 



For any a and b satisfying \a\, \b\ < 1 it is well known that 

\a N -b N \ < N\a-b\. (A24) 

Applying this inequality to our problem gives 



8(t,N) < N\(l-^-K(t/VN))-(l 



AN' 



< \pt 2 (K(t/s/N) - l)|/4 



(A25) 



For fixed t, the argument of K(t/s/N) vanishes as N in- 
creases, entailing pointwise convergence. Consider a point 
t, for any e > we can always find a u e such that \K(u) — 
1| < Ae/(vt 2 ) for all |u| < u e . The variable u t defines a 
N e = t 2 u 2 , such that \u\ = \t/y/N\ < u e whenever N > N £ 
and so 



S(t,N>N e ) < e. 



(A26) 



Hence we have shown Thm.[3] Furthermore, this convergence 
is uniform on finite intervals of t. 



Appendix B: Convergence of expectation values 

In the letter we state that convergence of \a„ t° a Gaussian 
function with covariance matrix r CToo entails convergence for 



the expectation value of operators, B, with absolutely inte- 
grable characteristic functions. Here we prove the relevant 
lemma, and also show that its conditions are meet for opera- 
tors of the form B = \x)(y\ where \x) and \y) are, up-to a 
Gaussian unitary, multimode Fock states. Note that, when a n 
are positive Hermitian operators, stronger conclusions can be 
reached by a similar means as shown in Ref. [ 1 1|. 

Lemma 3 (Convergence of expectation values) Consider a 
sequence of trace class m-mode operators a n and a limit- 
ing operator a^, such that (i) \Xa-„ oo( r )| < 1 for all r 
and (ii) For any compact region, R, \a n converges uniformly 
to x&ce- Consider also a trace class operator B satisfying 
J |xs(r)|rfr < oo and any e > 0. There exists an n e such 
that for all n e < n 



\tr(Ba n ) - tr(Ba Q 



< e. 



(Bl) 



For the sequence of characteristic functions considered in the 
letter, property (i) is a prerequisite and property (ii) was shown 
in Cor. ([T}. We state these properties again here for clarity. 
First we note that for any two m-mode trace class operators, 
B and A, the characteristic functions, \b and xa, satisfy 

J dr XB (r) X A(r) = (2n) m tr(BA). (B2) 

This entails that 

D n - (27r)" l |tr(B ( T„)-tr(B ( 7 co )| (B3) 

XB (r)[x,Ar) - xa x (r)}dr\. (B4) 



Next we observe that for any xb satisfying j \xb (r) \dv < oo, 
and for any e > we can find a function h with compact 
support such that 

J \ XB (r) ~ h(v)\dr < e. (B5) 

That is, the set of functions with compact support are dense 
in the set of absolutely integrable functions. Defining 7j(r) = 
Xs(r-) - h(r) and A„(r) = x<x„(r) - Xa^i^), we have 



D n = | [r](r) - h(r)]A n (r)dr\, (B6) 

< | J ??(r)A„(r)dr| + | J /i(r)A„(r)dr|, 

< sup(|A„(r)|) [\r,(r)\dr+\ f h(v)A n dr\, 

< 2e+| f h(r)A n dr\. 

The third line uses a Holder's inequality. The final line uses 
Eq. (B5i and |A n (r)| < 2. For sufficiently large n, the sec- 
ond term can also be made arbitrarily small because h has 
compact support and A„ vanishes uniformly on any compact 
set, as required by property (ii). This completes the proof of 
the lemma. 



Furthermore, any operator B — \x)(y\ will satisfy the con- 
ditions of the lemma when \x) and \y) are eigenvectors of a 
Gaussian density operator. That is, there exists a Gaussian 
unitary such that U\x) and U\y) are Fock states. Clearly B is 
trace class, and so we only need show absolute convergence of 
Xb- For a 1-mode element in the Fock basis \n) (m\ the char- 
acteristic function will be a product of a real polynomial in r 
of finite degree and a decaying Gaussian in |r|. Such functions 
are easily seen to absolutely converge, and so too shall a prod- 
uct of m such functions. A Gaussian displacement only adds a 
phase factor to the characteristic function, and non-displacing 
Gaussian unitaries are equivalent to a unitary change of vari- 
ables in the integral. Hence, Gaussian unitaries will not alter 
the value of the absolute integral. Indeed, most trace class 
operators are easily seen to have absolutely integrable charac- 
teristic functions. However, we do not currently know if this 
is generally true. 



Appendix C: Gaussian integrals and the Schur complement 

The product of two operators, e.g., pll has a characteris- 
tic function that equals an integral involving the characteris- 



tic functions for p and n as expressed by Eq. ( 15 1 of the let- 
ter. Furthermore, for Gaussian states this integral is solved 
by taking the Schur complement of an appropriate matrix (see 
Eq. ( [To} of the letter). These techniques are standard, but for 
the purpose of completeness we shall review them here. 

As a side remark, note that an alternative strategy to- 
wards formulating generalized quantum central limit the- 
orems would have been to consider objects of the form 
U 1 / 2 pU 1 / 2 /tr(np) instead of IIp/tr(IIp). While this con- 
stitutes an alternative in principle, the resulting Schur com- 
plement expressions become much more involved, motivating 
the path taken above. 

We begin by noting that for any trace class operator A, the 
formula for the characteristic function can be inverted such 
that 

A ex f dr XA (r)D(-r). (CI) 

Hence, the product of two operators satisfies 

ABcx f f drdr' X A(r)xB{r')D{-r)D{-r'). (C2) 



The Campbell-Baker-Haussdorf formula entails that 

D(-r)D(-r') = D(-(r + r')) cxp(-ir T Sr72), (C3) 
and so 

AB ex /drdr'xA(r)xs(r')£'(-(r + r'))exp(-ir T Er72). 

(C4) 
Changing variables to q = r + r' and r = r yields, 

AB cc drdqxA(r)xB(q-r).D(-q)exp(--r T I](q-r)). 

(C5) 



Noting that r T Er = for all r eliminates part of the exponent. 
Next we observe that this has the form of Eq. (|CT]) and so we 
can deduce that the characteristic function is 



XAs(q)oc / drxA(r)xB(q-r)cxp(--r Eq). (C6) 



Setting A = poo and B = U provides Eq. 15 of the letter. 



Until now we have not used any properties of A or B, but 
herein assume they are Gaussian with covariance matrices T A 
and T B . Hence 



XAs(q) oc / drexp 



(C7) 



x exp 



(q-r) 7 r B (q-r) 



-r J Eq 



4 2 

This can be expressed more compactly using a single matrix 

XAB(q) « y C xp(-(q,r) T Af(q,r)/4)dr, (C8) 

where M is the Block matrix 

M = x 2iI-r B Tllr B )- «*> 

Observing that rEq = — qEr allows explicit symmetrization, 
T B -T B - ffl 



M = 



(CIO) 



ffi - r B t a + r B 

such that now M = M T . The Schur complements arises by 
decomposing the matrix as M = NDN T , with 



D 



Sc(M) 

T A + Y L 



(Cll) 



;v = [ J-^ + ^ + r,)- j (C12) 



where Sc(M) denotes the Schur complement of M, 

Sc(M) = T B - (T B + iE)(T A + T B )-\T B - »E). (C13) 

This decomposition requires only that M is symmetric and 
that T A + T B is invertible. Returning to our Gaussian integral, 
a change of variables (q, r)N T H> (q, r) gives 

XAB(q) oc Jexp(-(q,r) T D(q,r)/4)dr, (C14) 

and so the q and r variables decouple, and the integral evalu- 
ates to some number and 

Xj4 i3(q) ex exp(-q T Sc(Af)q/4)dr. (C15) 

Hence, we have that T AB = Sc(M). By substituting in T A — 
r Poo and T B = Tn we obtain Eq. ([T6|) of the letter. 
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Appendix D: Notions of weak convergence 

Here we briefly discuss technicalities related to notions of 
weak convergence versus convergence in trace-norm. Interest- 
ingly, these different notions can have a physical implication. 
While in many experimentally relevant examples the condi- 
tions of Thm. Q of the letter are satisfied, one can concoct 
exotic instances for which strong convergence fails. For ex- 
ample, if we consider the 2-mode state vector 



hfo>«|0 1 0)+0.1|l,l>+6|2,2) 



(Dl) 



and apply the GP protocol, then Thm. (J2]i of the letter holds 
for a target Gaussian 



IV>C 



Ea/ior 



n, n) 



(D2) 



However, once correctly normalized, straightforward numer- 
ics show the state diverges with fidelity (t/j n \'4'oo) vanishing 
with n. This highlights the importance of proving more ro- 
bust convergence results. Stronger convergence has only been 
established for EP, where trivially tr(p ra n) = tr(/o oo n) since 
II = t. On a theoretical level, it is compelling to ask what the 
necessary and sufficient conditions are for strong convergence 
ofGG 



Appendix E: Iterative formula 



Here we go through the details of the iterative formula of 



Eq. ( 19 1. For brevity, we assume n is diagonal in the Fock 
basis, and so V = 11, though the more general case can be 
proven by the same method but more cluttered notation. For 
the (n + l)th round we have that the state satisfies 



x,y 

T n+1 ~ 



tr 2 [U(p n ® p n )W(IL®IL)] 
ix[U{p n ® p n )W(JL®U)] 
ti 2 [U(p n ®p n )W(IL®IV)} 



tr[p„IT] 2 
and so the expectation values satisfy 

xy ti[(® k a x k '; 1 )H® J a v / A )U(p n ®p n )W(U®U)} 

a " +1 ~ ' tr[ P „np 

= ti[U\® k a x k k ^ {® ] a v j 3 l )U{a n ® a n )}. 

Next we must find the effect of conjugating the unitary 
through the annihilation/creation operators. 



^Ec^i®^^®^-^®^^®^-^) 

where in the last line we have simply noted that all the bi- 
nomial factors and contributions of l/\/2 are non-negative 
and real numbers, and so sum to a non-negative number 



C*y. The exact value of these numbers are unimportant to 
us here. We substitute this expression for the conjugated anni- 
hilation/creation operators into that for a*^ . Evaluating the 
trace, we replace the operators on the right hand side with the 
variables of the form a" ,v and arrive at the desired formula. 



Appendix F: Trace norm convergence 

Here we give our own version of the well known proof of 
trace norm convergence of p n to p^ from point-wise conver- 
gence of the density matrix elements. First we note that for 
any e > there exists a projector P of finite rank that com- 
mutes with II, such that tr(Pp oc ) > 1 — e, It is useful to 
denote P — J^xex \ x )( x \ wnere X is a finite set of eigen- 
vectors of n. Furthermore, we can find an n e such that for all 
ii > n f we also have 



tr(Ppoo) - tr(Pp n ) < e. 
Combining this with tr(Ppoo) > 1 — e, gives 
tr(Pp n ) > 1 - 2e. 



(Fl) 



(F2) 



Next we observe that the trace norm can be bounded as fol- 
lows, 

\\Pn - Poo||l 
= \\(P(Pn - Poo)P + (Pn - PPnP) ~ (Poo - P PooP)\\l 
< WPiPn-p^PWl + WPn-PpnPWl 

+ \\Poc - PpooP\\l 

(F3) 
Let us consider the first term, which satisfies 

\\P(Pn~Poo)P\\l< Y, \(y\Pn\x) ~ (y\ P oo\x}\. (F4) 
x,y£X 

By our previously established results, for every pair \x),\y) e 
X, there exists an n eyX ^ y such that for all n > n e , x , y , we have 

\(y\p n \x) - (y\ Poo \x)\ < e/\X\ 2 (F5) 

Hence, for all n > n* where n* = iJia,x{n eXt y, n e } we have 

\\P(p n ~Poo)P\\i<e. (F6) 

We next consider the third term of Eq. (|F3|l, noting that 



\\Poo - PpooP\\i < 2\\P Poo (l-P)\\ 1 , (F7) 

+||(l-P) Poo (l-P)||i. 

The later term is simply ||(1 - P)p 00 {l - P)||i = tr((l - 
P)Poo) < e. Since p^ is Hermitian and positive we deduce 



WPfaO—P)^ < y/\\P Poo P\\ 1 \\(l-P)p 00 (l-P)\\ 1 , 

which easily shown to be an equality when p^ is pure, 
with the inequality following for mixed states. Noting 
\\PpooP\h < land again ||(1 -P) Poo (l-P)|| 1 < e gives 

||Ppoo(l-P)||i <V~C 
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and so 

\\Poo~PpooPWi < 2V~e + e. (F8) 

By a similar argument we can also show that for all n > n e , 

\\pn-Pp n P\\i<V2~e + 2e. (F9) 

Putting these pieces together entails that for all n > n* , 

||pn-A»||i<4e+(2 + V2)V£. (F10) 



Hence, for any desired trace norm accuracy, and hence also 
fidelity, there always exists a number of iterations, n*, above 
which the protocol achieves this accuracy. 
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